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Abstract 

In this paper we treat the so-called dynamical Casimir effect in a cavity with a 
two-level atom and give an analytic approximate solution under the general setting. 

The aim of the paper is to show another approach based on Mathematical Physics 

to the paper [arXiv : 1112.0523 (quant-ph)] by A. V. Dodonov and V. V. Dodonov. 

We believe that our method is simple and beautiful. 
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1 Introduction 



In this paper we revisit the so-called dynamical Casimir effect (DCE). This means the photon 
generation from vacuum due to the motion (change) of neutral boundaries, which corresponds 
to a kind of (quantum) fluctuation of the electro-magnetic field. 

This phenomenon is a typical example of interaction between microscopic and macro- 
scopic levels and is very fascinating from the point of view of not only Physics but also 
Mathematical Physics. See for example pQ and its references. 

Then, how do we detect the photons generated ? This is also an important problem. 
Recently, Dodonov and Dodonov in [2] treated this problem by use of a two-level system of 
an atom in a cavity. They called it "the cavity dynamical Casimir effect in the presence of 
a two-level atom", and gave an effective Hamiltonian by simplifying the Hamiltonian given 
by Law [3] and adding a two-level system of an atom to it, and constructed an approximate 
analytical solution. 

In this paper we treat this problem once more and present another approach to construct 
an analytic approximate solution under the general setting. We believe that our method is 
clearer than that of [2]. 



First of all let us make a brief review of [2J within our necessity. In the following we set 
h = 1 for simplicity. 

(i) Cavity DCE ([3], [I]). For this we take the simplest Hamiltonian which is the special 
case of Law [3] (namely, e(x,t) = e(t)) 



where u(t) is a periodic function depending on the cavity form and x{t) is given by 



2 Model 



H DCE = u{t)a)a + ix{t){{a)) 2 




(1) 



x(t) 



1 dw(t) 1 d 



log|w(t)| 



4w(t) dt Adt 
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, and a and cr are the cavity annihilation and creation operators respectively. Therefore, 
physics that we treat in the following is two-photon generation from the vacuum state. 



(ii) Detection ([1], [5]). For this we take a two-level system of an atom inserted in the cavity 
and the Rabi Hamiltonian as interaction 

H D = ^03 <g> 1 + g(a+ + 0_) <g> (a + a)) 

= -^cr 3 <g> 1 + #(0+ <g> a + a + (g) a 1 " + a_ ® a + a_ ® a 1 ") (2) 

where 1 is the identity operator on the Fock space generated by {a, at,N = a^a} and 

i o \ / o i \ \ 1 

0-1/ V 00 / V 10 / V 01 

However, to solve this system is very hard. In fact, our aim is to obtain some analytic 
approximate solution. Therefore, as a rule we replace the Rabi Hamiltonian with the Jaynes- 
Cummings one 

H D = ^(T 3 ®l+g((T + ®a + <T-®a i ). (3) 
Here, let us note the well-known su(2) relations 

[(1/2)03,0+] = 0+, [(l/2)0 3 ,0_] = [0+,0_] = 2 X (1/2)03- 

(iii) Total System. By adding fl3]) to (QQ) we treat the (effective) Hamiltonian like 

H = H[t) = Hqce + Hd 

= w(t)l 2 ®JV + ix(t)l2® {(a t ) 2 -a 2 } + ^0 3 ®l + 5'(0+(8)a + 0-«)a t ). (4) 

Our aim is to solve the Schrodinger equation 

• ^|*(t)) = J ff|*(t)) = J ff(t)|*(t)> (5) 

under the general setting. 

If H is time-independent, then the general (formal) solution is given by 

|tf(t)) = e - itH \v(o)). 



To calculate e~ ttH exactly is another problem (which is in general very hard). However, in 
our case H is time-dependent, so solving (J5J) becomes increasingly difficult, so we must give 
a further approximation to the Hamiltonian. 

(iv) Interaction Picture. From here let us change to the interaction picture. Namely, for 
V = V{t) we set 

m)) = yt|*(t)) Mt)) = vmt))- (6) 

Then it is easy to see that the equation (jSJ) can be changed to 

i±m)) = (vBv- t ^±v)m)). ( 7) 

This is in general called the interaction picture. 

In the following we impose some restriction on our model. Namely, we take u(t) as in [2] 

uj{t) = cu (l + esm(r)t)) 

where u , e and rj are some constants. We assume that e is small enough comparing to u . 
Then u(t) ~ u and 

ericos(rit) en . , 

x (t) = —-!- ; 7 « -/ cos (rft) 

AW 4(l + esin(^)) 4 KIJ 

, and here we take V as 

V = V{t) = e"^ 3 <g> = 



Then some calculation by use of ([HD gives 

Hit) = V*HV-iV*^V 

= (wb - |) l a ® JV + h ® K^) 2 - a 2 + e 2 "V) 2 - e~ 2 ^a 2 } 
+ g> 1 + £(a+ g> a + a_ g> a T ) 

and we also use the rotating wave approximation e 2ir)t ~ if 77 is large enough. As a result 
we finally obtain the time- independent Hamiltonian (see [2]) like 

H=(u - J) 1 2 ® iV + i^l 2 <g> {(a f ) 2 - a 2 } + ^ n 03 ® 1 + 0(0+ ® a + <r_ <8> a f ). (9) 
V 2 ' o 2 
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(v) Our Target. The aim of this paper is to solve the equation 



dt 

under the general setting and the formal solution is given by 



4l*(*)> = #!*(*)> ( 10 ) 



|$(t)) = e- tU1 |$(0)>. (11) 

Therefore, what we must do in the following is to calculate the term e~ ltH explicitly. 
Before closing this section we must present an important problem. 

Problem We used approximation (a kind of rotating wave approximation) two times. Make 
an adaptive range (or region) of the model clear. 

3 Approximate Solution 

In this section we try to calculate the term e~ ltH in (jlll) . For simplicity we set 



-UH = X + Y 

where 

X = -itA, A = (lo - |) h ® N + i^-h ® {(at) 2 - a 2 } , 
— r)/2 

Y = -itB, B = Y~ 0, i ® 1 + 9( a + ® a + or - ® flt )- ( 12 ) 

For our purpose the Zassenhaus formula is useful : 
Zassenhaus Formula We have an expansion 

e X+Y = . . . e -\{n[X,Y],Y]+[[xy\,X}} e \[X,Y] e Y e X _ 

Note that the formula is a bit different from that of [6] . In this paper we use 

e X +Y ^ e \[X,Y] e Y e X = e -* [A ,B] e -itB e -itA_ (13) 

Therefore, let us calculate each term in the following. 



[I] First, we calculate e ltA . If we set 

a = 2(wo-![), = \en (14) 

then A in (Tl2l) becomes 



A = |l 2 ® AT + i|l a ® {(at) 2 - a 2 } = 1 2 ® jc^iV + z/3 Q( a t) 2 - l -a 2 
Here, we use a well-known Lie algebraic method, see for example [TJ. Namely, if we set 

#+=^ f ) 2 , *- = ^ 2 , K * = l( N+ l) (15) 
it is easy to see both K\ = K_, K\ = K 3 and the su(l, 1) relations 

[K 3 , K + ] = K+, [K 3 , K-) = -K-, [K+, K-) = -2K 3 (16) 

by use of the relation [a, a'] = 1. Then, A in ( |12p can be written as 

A = l 2 ®^-~a + aK 3 + ip(K + -K_) 
-\a + aK 3 + i(5 (K + - KJ) 

-\a + aK 3 + i(3 (K + - KJ) 

and from this we have 

( it a -it{aK 3 +iP(K + -K-)} 
e I±a e -it{aK 3 +if3(K+-K-)} 

Therefore, we have only to calculate the term 

JJ{+\ = e -it{aK 3 +ip{K+-K-)} _ e -itaK 3 +t0{K + -K-) 

For the purpose we want to look for the following disentangling form 

U(t) = e fm +e 9(t)Ks e h(t)K- (17) 

with functions f(t), g(t), h(t) (/(0) = g(0) = h(0) = 0)Q. 



this is a standard method 



The result is as follows. 



f(t) 



4^ 



£ sinh I tJ-^ + fi 2 



cosh [tJ-^+P 2 + 



— sinhlW-f + /3 2 



g(t) = -21og(cosh(^/- — + /3 2 ) + 



a' 



+ /3 2 



: sinh t 



+ /? 2 



p 2 sinh + ^ 



COsh UJ-f + /3 2 + 



1 sinh y-f + /3 2 



-/(*)• 



The derivation is analogous to that of [81. See also the appendix. 



[II] Second, we calculate e The result is well-known and is given by 



-itB 



exp < —it 



2 



ga 



( 



\ 



cost-^/ 'if + g 2 — 



— ig — 



r+9 2 a \ 



■ sin t.ftp + 



cost^+f^ y 



where we have set 



for simplicity. 



5 = n - -, if=—+g 2 N 
2 Y 4 



[III] Third, we calculate e~~^ A,B \ From (fT2l) simple calculation gives 



[A, 5] = ^^(-<r+ ®a + o-_ tgJa 1 ") -ig(3(a + ® a ] + a_ <g> a) 

-fa-i^at 
§a f - z/3a 



-L> 



Dt 



where we have set 



for simplicity. Note that 



D = —a + iBa\ = —a) — iBa 
2 ' ' 2 ' 



The result is 



e 2 



[A,B] 



cos (^v/^Dt) S i n (g±y/m) D 



[IV] As a result, our approximate solution to the equation <^ is given by 
under any initial value |\l/(0)). If we write 



(t) = e 



then each component is given by 



i 1 



-\[A,B] e -itB e -itA 



U21 U22 



iS sint-^/ f + g 2 



U u e> { cos ( g-^/DDl) ( cos*vV + 9 2 ~ 2 - 



_ it 

C/12 = e- ( 



^sinf/V^tV^^atj 

\ 2 ; vVT^ 



e f(t)K +e g(t)K 3e h(t)K. 



-ig cos 



1 



t/21 = e 4 



sin g—V DD^ D costJf + _ 
sin J g^VDtZ^ D f ( cos ty/ip + g 2 



idsmt^f\ I e f{t)K +e g(t)K 3e 



id sint^f + g" 2 

Vdw~'"V 2 7" y"""* r ' a 2 ^/f~T7~ 



U22 = e* a < ig 



1 ft 2 r^—\ + sini v / ^+ 9 2 
^= sin g-VDW D ] yy y a + 

vww V 2 ; VfT¥ 



t 2 1 — 7 — \ / ^_ i5sini,/<2j 
f ,, s , g-VDW j \wety/ip + - 



e f(t)K +e g(t)K 3e h(t)K. 



This is our main result. 
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4 Closing Remarks 

In this paper we treated the model by A. V. Dodonov and V. V. Dodonov and presented 
another approach based on Mathematical Physics and obtained some analytic approximate 
solution. Details with applications and developments will be published separately. 

We have neglected Dissipation (for example, the Cavity loss) in this paper, which is not 
realistic. If we take dissipation into consideration the model will become (very) compli- 
cated. For example, see our papers [H], PH] and [TT], [12]. Then, to obtain an analytic 
approximate solution will become increasingly difficult. 

Further study and new ideas are needed. 



Appendix 



A note is added. In the text we must calculate 



U(t) = e 



-it{aK 3 +i/3(K+-K-)} _ -itaKa+tP(K+-K-) 



From this we have the differential equation 

d 



dt 



U(t) = {-taK 3 + p(K+ - K.)} U(t). 



On the other hand, under the ansatz (fT7j) (: U(t) = e^ t ' K+ e 9 ^ t ' K3 e h ^ t ' K ^) some calculation 
gives 

jU(t) = {(f-gf + he-vf) K + +(g- 2he^f) K 3 + (he-*) K_}u{t) 
where we have used / = %, etc for simplicity. Therefore, by comparing two equations we 



obtain 



f-gf + he-9f 2 = f3 
g — 2he~ 9 f = —ia 
he~9 = -f3 



f-9f-Pf 2 = P 
g + 2(3f = -ia 
'he' 9 = -f3 



g + 2/3f = -ia 
h = ~(3e g . 
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The equation 

f + iaf + l3f = p 

is a (famous) Riccati equation of general type. If we can solve the equation we have the 
solutions like f(t) =>- g(t) =>- h(t). See (fl8|) as these solutions. 
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